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. Abstract 

We consider the asymmetric six-vertex model, i.e. the symmetric six-vertex model in an external 



field with both horizontal and vertical components, and the relevant asymmetric XXZ chain. The 



a 

q | model is widely used to describe the equilibrium shape of a crystal. By means of the Bethe Ansatz 



solution we determine the exact free energy singularity, as function of both components of the 
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experimentally), as the antiferroelectric ordered phase is reached from the incommensurate phase 
normally to the phase boundary, and we determine a new singularity along the tangential direction. 
Both singularities describe the rounding off of the crystal near a facet. The hole excitations of 
the spin chain at this point on the phase boundary show dispersion relations with the striking form 
AE ~ (AP)a at small momenta, leading to a finite size scaling AE ~ iV~2 for the low-lying excited 
states, where iV is the size of the chain. We conjecture that a Pokrovskii-Talapov phase transition 
is replaced at this point by a transition with diverging correlation length, but not classified in terms 
of conformal field theory. 

PACS numbers: 05.70.Jk, 64.60.-i, 64.60.Fr, 75.10.Jm 
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1 Introduction 



After the pioneering paper of Yang, Yang and Sutherland [T|] the asymmetric six- vertex model, 
i.e. the symmetric six-vertex model in a field, was recently rediscovered because of its con- 
nection to a number of physically interesting problems, first among them the determination 
of the shape of a crystal at equilibrium with its vapor phase ||, || . This can be achieved by 
mapping the asymmetric 6-vertex onto, say, the (001) facet of a bcc crystal under the condi- 
tion that no overhangs or voids are allowed (see e.g and |4[] for details on the mapping). 
Excitations in the vertex model correpond to small tilts away from the (001) facet @] and 
it can be shown || that the free energy as function of the two components of the field gives 
exactly the equilibrium shape of the crystal. 

In its own, the asymmetric six-vertex model provides an interesting 2-dimensional system 
of interacting dipoles in an external field with horizontal and vertical components (h,v). 
Fluctuating two-valued variables (dipoles) are attached to the links of a two-dimensional 
square lattice, and the model is defined by assigning a set of Boltzmann weights (equivalently, 
interaction energies) to each allowed vertex configuration (see fig. 1). The transfer matrix 
can be diagonalized exactly by the Bethe-Ansatz in its coordinate or algebraic version || fzj. 
The phase diagram and the nature of the phase transitions are well understood when h = 
v = (symmetric six- vertex), or when h = and v ^ ||, |J. If h, v ^ 0, some general 
features of the phase diagram have been described in |I] and the details of the calculation 
spelled out in ||, but a few questions have remained unanswered. It is known that, in the 
antiferroelectric regime, with which we will be concerned in this paper, the free energy f(h, v) 
remains constant as function of the field ('flat phase') in a bounded region of the (h,v) plane 
containing h = v = 0. This corresponds to the flat (001) crystal plane. Beyond this region, 
bounded by a curve T, the field is sufficiently strong to destroy the antiferroelectric order 
of the system, but not strong enough to impose ferroelectric order, and an incommensurate 
phase appears where the polarization (zero in the flat phase) changes continuously with the 
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field. Here the spectrum of the transfer matrix is gapless with finite size corrections typical 
of the gaussian model |n|. The curve T has been investigated only at a few points and 
the nature of the phase transition was found to be of a Pokrovskii-Talapov (PT) type |TT| . 
In this paper we determine the exact free energy singularity when approaching T from the 
incommensurate phase, which determines the curvature of the crystal near the (001) facet. 
We find two special points (h c , v c ) (by symmetry under arrow-reversal the same holds for 
(—h c , —Vc), also on T), where 

3 

f(h c + 5h,v c ) = f(h c ,v c )-const(5h)z 

f(h c ,v c + 5v) = f(h c ,v c ) - const\5v\ 3 (1.1) 

These are points on T where the tangent to T is parallel to the v— axis. The exponent 3/2 
measures the rounding off of the edges of the (001) facet. 

Even though our calculation has been carried out only at these points of T, the technique 
we present should work in general, and our result, which generalizes Lieb and Wu's method and 



complements the finite-size techniques of Kim [10], strengthens the long held belief that the 



exponent 3/2 should govern the free energy singularity at every point of the phase boundary 
HH, p. This exponent has been measured |T^| in some experiments with Pb crystals some years 



ago. Our results however show that along some tangential direction the exponent 3 should 
dominate, and this fact should be observable experimentally. 

However, something more can be said about the nature of the phase transition along T. 

The method of mapping a 2d statistical system into a Id quantum spin chain has been 



fruitful and widely used in the past |13|. We pursue it here, regardless of the fact that the 



relevant spin chain which turns out to be the asymmetric XX Z spin chain in a vertical field 



V, is not hermitian [Q. We find that the flat phase corresponds to a region in the (h,V) 
plane where the ground state energy does not depend on the fields and where excitations are 
massive. Along the transition line, analogue of T, the excitations become massless, but the 
point (h c , V — 0) (with its symmetric (— h c , V = 0)) is singled out by the fact that dispersion 
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relations obey the striking law, at small momenta 



AE ~ (AP) 



1/2 



(1.2) 



and finite-size corrections for low-lying excitations scale like 



AE ~ N~ 1/2 



(1.3) 



if iV is the length of the chain. At this point, analogue of (h c , v c ) for the statistical model, 
the vanishing of the mass gap exhibits an exponent 1/2 which does not appear at any other 
point of the transition line in the (h, V) plane. In section 5 we propose an explanation of 
these results, arguing that at the point (h c , v c ) the transition occurs with a divergence of the 
correlation length in the correlator between two vertical arrows, while everywhere else the 
transition is induced by level- crossing (Pokrovskii-Talapov) which prevents the divergence of 
the same correlation length. 

The paper is divided in 5 sections. In section 2 we give definitions and summarize pre- 
viously known results. In section 3 hole excitations and the spectrum of the spin chain are 
studied and in section 4 the method of Lieb and Wu is suitably extended to determine the 
free energy singularity when both h and v are nonzero. Section 5 contains an interpretation 
of the results. 



The model is a natural generalization of the well-known symmetric six-vertex model. Arrows 



2 Definitions 



are placed on the edges of an N x M square lattice and Boltzmann 




assigned to the vertices (see Fig. 1) so that the row-to-row transfer matrices 





form a commuting family 



[T(«),T(«')] = 



4 



for any two values u, vl of the spectral parameter [|15|]. The associated (integrable) spin chain 



cosh 7 ^ __ ; ^ 2; , , _, lh 



e (T, (J; 



(7 \ 



[2.2) 



is obtained from ( |2.1|) by taking the so-called extremely anisotropic limit (u — > 0) 



(«+fc)E?=i°?7Y^ ft = -F-^logfe^Ef^J) - s inh7^-logT(w)| u= o 



V breaks the Z 2 symmetry of spin reversal while h breaks parity invariance (see App. A for 
a complete discussion of symmetries). 

By means of the Bethe-Ansatz, eigenvalues of ( |2.1| ) and ( |2.2| ) are found from the solution 
of a set of coupled equations 



sinh(j + ^) l " = +12hN " sinh( 7 + |( afc - a/ )) 

sinh(f-^)J " l J i i 1 sinh(7-|K-aO) 



A; = 1,2, 



(2.3) 



and given respectively by 
A(u) 



e ^(AT-2n) e / l Ar 



sinh(7 — w)]-^ " sinh(| + « — 



sinh7 



2 



J'=l 



sinh(| - u + ^) 



+e ^(V-2n) e -hV 



sinhwl^ " sinh(^ + u-^ 
II 



sinh7 



4 «i 

2 
2 



=1 sinh(| - w + ^-) 



= A R {u)+A L {u) 

£ = iVcosh7-y , 7 V(N-2n) 

r^, cosh 7 — cos 

n 

= iycosh7 + ^e(a fc )-\/(A^-2n) 

k=l 



(2.4) 



(2.5) 



Here n stands for the number of reversed spins (arrows) with respect to the reference ferro- 



|). It is a conserved quantity since S z = \ Y^=i a j commutes with Tt 
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magnetic state \j\ 
and T{u). 

Beside their energy, given by fl2.5p, the momentum can also be computed. T 1 (0) yields 
the right-shift operator S = e _iP 



S = |«i, a 2 , • • • ,«m) = |aM,ai,«2, • • ',ac M -i) 



(2.6) 
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and from ( |2.1|) and ( |2.4| ) one gets 



-iP 



with 



p°(a) 



-2 hn sinh(l + 
j=1 sinhQ - Si) 

sinh(§ + f 



P = -2m/j - ^2p°i 



(2.7) 



-i In 



.sinh(§ - f ) 

Unlike most integrable spin chains studied before, (|2.2| ) is not hermitian for h ^ 0, even 
though the statistical model has physically sensible positive Boltzmann weights. The question 
arises whether (|2.1| ) and (|2.2| ) have a complete set of eigenvectors. By numerically diagonal- 
izing the transfer matrix on small chains it appears that, even though a few eigenvalues are 
degenerate in some charge sectors, the eigenvectors are linearly independent. We will assume 
that a complete set of eigenvectors exists and it is given by the Bethe Ansatz. 

We first summarize, for the sake of completeness, some already known facts that have 
been published elsewhere, beginning with the spin chain because, although the Bethe Ansatz 
equations are identical, the form of the energy contribution from each rapidity a^ is simpler. 
Taking the logarithm of (|2.3|) in the usual way we get 



1 n 2n 
p°(a k ) - — ®( a k - <*i) + 2ih = -j^h 



k = 1,2, 



n 



[2.8) 



where Ik is half-odd (integer) if n is even (odd) and 
-sinh( 7 + f ) 



9(a) 



-i In 



_sinh(7 - y) 

We define p°(0) = O(0) = and cuts are chosen to run from 27 to zoo and from — 27 to 
—ioo for p°(a); from 2i / y to too and from —2^7 to — ioo for 0(a). Notice from ( |2.3| ) that 
Re(a) G [— 7r,7r] so Re(ak — on) G [— 2n, 27r]. The cuts are chosen so that 0(a) is analytical 
in —27 < Im{oc) < 27 and real, monotonically increasing when a G [— 27T, 2ir]. With these 
conventions, the ground state at h = in each sector of fixed S z = y — n corresponds to a 
sequence of n consecutive numbers {Ik} in (|2.8j), from — to symmetric around [ jlBf . 
The rapidities {aj} are real and distributed symmetrically around a = too. As h 7^ 0, the 



rapidities move into the complex plane along a curve C. In a standard way |9[], one gets in 
the thermodynamic limit from 



PV) - 77- / ^@(« - P)R(P) + 2iA = 2ttx - ^— ^ < x < — (2.9) 
2ti Jc 4 4 

x is the real parameter of the curve, y is the polarization defined through 

y = lim = lim (1 - =-7) (2.10) 

and the rapidity density R{a{) = linijv^oo j^t^z^s * s determined by solving the integral 
equation 

£(a) / d[3K(a - p)R{8) = R(a) (2.11) 

2tt Jc 



where 



£(«) = #(a) - rf0(a) 



da ' da 
The energy and the polarization are thus given by 
E 



lim 

N^oo N 



- cosh 7 + — / dae{a)R{a) - Vy (2.12) 
2ir Jc 

l ~ir = s/ c rfaJi(a) (2 - 13) 

Some preliminary information about the shape of C and its location in the complex plane can 
be obtained by solving ( |2.3[ ) numerically. We take as initial solution that composed of real 
roots and corresponding to the ground state, at fixed S z , for h = 0. By the Perron-Frobenius 
theorem |16], [17|], the relevant eigenstate remains the ground-state at fixed S z even when 
h 7^ 0, and it is real. It appears that for h > (h < 0) the rapidities move into the lower 
(upper) half-plane as shown in fig. 2. 

The curve C is invariant under a — > —a*, which is to be expected, being this transforma- 
tion a symmetry of (|2.3| ), so we will set A = —a + ib and S = a + z& to be the endpoints of 
the curve. Note that this property makes E real, as it should. 



Strictly speaking R(a) is defined on C only, but ( |2.11| ) can be used to define it outside of 
C. If C is contained in the strip —7 < Im(a) < 7, R(a) is analytic in —7 < Im(a) < 7, but 
it inherits the poles of at ±27. Let us consider the curve for which a = tt. In this case, 



since R(a) is 27r-periodic, (|2.11|) can be solved straightforwardly by Fourier transform. The 
solution is 



fl(«) = E- 



— 7 < Im(a) < 7 



(2.14) 



2 cosh •yn 

Here and in the following, sums are understood to run from —00 to 00 unless otherwise stated. 
Beyond this strip ( |2.14j ) can be expressed using elliptic functions. Introducing the complete 
elliptic integrals of the first kind I (/') of modulus k (k r ), with k 2 + k' 2 = 1 |I8 |, related to 7 
by 



I{k) 



1 

71 



the solution of ( [2.1 1|) in a wider domain reads [19|, [18 



R(a) 



I(k) , fl(k)a 
-dn' 



7T 



\k 



71 



(2.15) 



Notice the presence of a pole at a = ±27, inherited from ^(a), which prevents the convergence 
of ( p. 14 ) beyond the smaller domain. The energy remains constant at its value for h = 



e 



hm — 

N^oo N 



cosh 7 — 2 sinh 7 V] 



-7|n| 



2 cosh jn 



and from Q2.13 ) y — 0. In fact the solution considered here has 



n 



iV 



rapidities (S z = 0) and, 



as it will be shown in the next section, E is the ground state energy for h and V sufficiently 
close to 0. 

As to the precise position of the curve, one has to revert to ( |2.9| ). Since 9 (a + 2n) = 
Q(a) + 2tt, we use the expansion 

g— ina— 2-y\n\ 



9(a) = a + i 



11 



27 < Im(a) < 27 
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and 

ina— y\n\ 

p°( a ) = a + iY- (2.16) 

and we introduce 

, N a _ e~ ina /I(k)a ,\ . 

p(a) = - + iY : =am(-^— ;k) 2.17 

yK ' 2 ^ 2ncosh 7 n If'/ v 7 

The series in ( |2.16 ) and ( |2.17 ) are certainly convergent when —7 < Im(a) < 7 but they 



converge also at a = ±n ± ij, because of the alternating sign. Eq. fl2.9p reduces to 

ib e nb 1 1 

p(a) + - +iY(-) n + 2ih = 2nx - - < x < - 

FK 1 2 2ncosh 7 n 4 ~ ~ 4 

Specialization to the endpoints permits relating the value of h to b || [| 

h{b) = -*-£{-)» Shl \ nb (2.18) 
v ' 2 ^ ' ncoshn 7 v ' 

so that the final equation of the curve is 

p(a) + ih = 2ttx — — < x < — 

Notice that the points on the curve are characterized by 

Im(p(a)) + h = (2.19) 

We set h c = h(b = —7). One might suspect that, when h > h c , the endpoints would remain 
at a = tt but with b < — 7 (or b > 7 if h < —h c ). If C does not cross the point — 27 where 
£(a) has a pole, it can always be deformed to the real axis in ( [2.11| ) 



R(a) + T duK(a - u)R(u) = £(a) (2.20) 
2tt J-tt 

so that the solution is still given by -^dnf ^-a; k), but the expansion fl2.14|) is no longer 



useful. To find the h(b) relation, we close C in (j2.9|) to the real axis, and take a = A 

p°(A) - r du&(A - u)R(u) + [ d(3R((3) + 2ih = - J (2.21) 

2lt J-TT J-TT 2 
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The integral of R(a) is obtained by integrating both sides of ( |2.20| ), and from ( p.21| ) we 
conclude 



2k( b ) = ^Sln ^T " j + 2U-r e '^T hnb (2.22) 



2 1 2 

which reduces to (|2.18| ) when |6| < 7. Eq. (|2.22[) though cannot give the right dependence 



h{b) at |&| > 7, because h{b) decreases when b < —7 and increases for b > 7, going back to the 
range of values it had as b G [—7,7]. Clearly the initial assumption a = it cannot be correct. 
To gain more insight we resort as usual to the numerical solution of ( |2.3| ) which shows that, 
as h > h c , the curve with y = has endpoints at (see table 2. Note that all extrapolations 
presented in the tables have been done using data up to 80 sites. We give however only 
the first values, up to iV = 40, for they already show clearly that the values are converging 
towards a limit. All tables have been calculated for cosh 7 = 21, where exp(2/i c ) = 10.51787). 

b = —7 a < 71 

This (new) result will be used in the calculation of the free energy singularity in Section 4. 

Turning next to the statistical model, the largest eigenvalue of the transfer matrix yields 
the free energy per site (we drop here the inessential factor (3) 

f , , s r lnAo(u,7>M) 
j(u, 7, n, v) = — hm — 

whose value, as iV — > 00, is dominated by the largest of the two limits 

1 a / \ j-. / r x 71 sinh(7 — u) 

hm -A r m = F R (u, 7, h,y)+vy = h + ln — h vy 

n^co JM sinn 7 



1 

Yn Jc 



daR(a)fn(a; u) 



1 sinh(w) 
lim ^tAl(m) = Fl(u, 7, h,y) + vy = -h + la—r^r h vy 

Af^oo jv smh 7 



[ daR(a)f L (a;u) 
2n Jc 
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where we have defined 



f R (a; u) = In - 



sinh(^ + u — ^) 



sinh(| — u + 

sinh(^+^-^) 

/ L (a;u) = In — 

U ^ sinh(|- M + ^) 

If Fr, Fl are known, and we call F the dominant one, the equilibrium value of y and the free 

energy are determined by the minimum condition 

f(u,<y,h,v) = min{-F(u,<y,h,y)-vy} (2.23) 

-1<2/<1 

When — h c < h < h c and for small enough values of v, the state defined by C (a = it; —7 < 

b < 7) also yields the largest eigenvalue of the transfer matrix. The free energy [] 

00 g— 2771 

f(u, 7, h, v ) = — 2 sinh(nw) sinhn(7 — u) (2.24) 

^ n cosh 

is constant in a whole region of the (h,v) plane bounded by a curve T ('flat' phase). The 
parametric equation (h(b),v(b)) of T is given by fl2.20|) and by 
dF 

h fixed, y=0 

which can be explicitly computed 



2,(6) = 7 - | 7 - 2 M + 6| + 2 g tr*^ ~ b ~ 2 " - "Dl - 7 < 6 < 7 (2.25) 

^ n cosh 717 

The other half of the curve T (see fig. 3) can be recovered from the symmetry f(—h, —v) = 

f(h,v) (see Appendix A). Most of these results have been obtained elsewhere and we have 

presented them here only for the sake of completeness. It should be pointed out that the fact 

that the ground state energy does not depend on h is simply a consequence of the analogous 

property of the free energy. 

lr The only part in which our analysis differs from [p| is that Ar is exponentially larger for d < 7 — 2u and 



Al for d > 7 — 2u where id is the point in which C crosses the imaginary axis. A comparison of (2.24) with 
[P| should take into account the different normalization of the Boltzmann weights. 
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3 Hole excitations of the spin chain 

To understand the nature of the phase transition along T we turn to the calculation of the 
excitation energies, and we set V = 0, since the role of V is simply to shift the spectra at 
S z 7^ 0. As proven in appendix A it is sufficient to consider h > 0. 

A complete treatment of the spectrum should rely on the classification of all possible 
solutions of (|2.3|) . This is usually done in the framework of the string hypothesis, according 
to which complex rapidities (at h = 0) have an imaginary part which tends to well defined 
values in the thermodynamic limit. Exceptional solutions other than strings, but that still 
appear in complex conjugate pairs, can be handled in a similar way [^(J. Yet, from the 
numerical analysis of ( |2.3| ), it appears that strings do not survive at moderately strong values 
of h. Therefore we shall limit our calculation to the so-called hole excitations, that is holes 
in the ground state distribution of rapidities, occurring in sectors with S z > (n < uq = 
We introduce the counting function 

n°(ry) 1 n h 

z (<*> W) = - ^ E ©(« - «,o + <- (3-i) 

so that (|2.8|) is rewritten 
Z(a k ) = | 

Set vac = number of vacancies available for the quantum numbers {Ik}- With the usual 
hypothesis that Z(a) be monotonically increasing, we have 

AZ = Z(a)\ Re ( a ) =7T - Z(a)\ Re (a)=~7T = -j-j- (3.2) 

On the other hand, from ( [3.1|) , we find AZ — 1 — For the ground state n = no = y and 
so vac = no, i-e. the available vacancies are all filled. For n rapidities, n = n — r where 
r = 1, 2, • • • we have 

AZ = 1 — — vac = n + r 

so that no + r vacancies are partially filled with (no — r) i^'s, leaving Nh = 2r holes. 
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We will resort to the 'backflow method' 0]| J2I] in dealing with (|2j), ([U]) and (£7|) in 
the limit A" — > oo. The calculation differs slightly for the 2 cases r = even or odd, but the 
results are identical and we will present the case r = even only. If this is the case, then n = uq 
(mod2) and the quantum numbers of the excited state have the same oddness of the 
quantum numbers {/°} of the ground state. We assume that the r additional vacancies for 
{Ik} are placed | to the left and | to the right of the sequence (see fig. 4). 



We call {Z?,- 1 ^} the | additional rapidities at the left edge, {/3¥'} the | additional rapidities 



(2)- 



at the right edge and {a, } the N h holes. Then, for the ground state 



Zn(a) 



Z (a° k ) 



p°(q) 
2tt 

N 



2kN& 



a: 



a]) + 1 



7T 



(3.3) 
(3.4) 



and for the excited state, adding and subtracting the holes 



Z(a) 



p°(a) 
2tt 

1 



"0 



— — V 0(a - a,) + — — y Q(a - af h 
2itN o^rM ^ K J 



3=1 



2kN& 



Z(a k ) 



2nN 

h 

N 



E 



3=1 



71 



Z(a 



j{h) 

_i 

N 



(3.5) 



(3.6) 



As N -»• oo {f3f ] } -> A and -> 5. Subtracting 



(2)- 



from ( |3.6| ) and retaining terms of 



order j, (this is a standard Bethe Ansatz calculation) |2l] we find that 



satisfies 



— a, 
hm — „ tt 



j(a) + — J c d(3K(a-(3)j(P) 



1 



Z7r i=i 



+ 



r 

47T 



G(a - B) + G(a - A) 



(3.7) 



with 



A£ = J c dae\a)3(a)-ye(af ] ) + r -{e(A) + e{B) 
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AP 



ihN h 



N h 



c 



da£(a)j(a) + J^P 



or Ah), 



[a 



j 



3=1 



p°(A)+p°(B) 



Eq. ( p.7|) defines the analytical properties of j(a) in the complex plane. Since for < h < h c 
the curve C is contained in —7 < Im{a) < 0, j(a) is certainly analytic in —27 < Im(a) < 7 
and the curve can be closed to the real axis. Noticing, from (pO] ), that j(a + 2n) —j(a) = 
we get, with u e [— it, n] 



j{u) + 



AE 



AP 



2n 



dvK{u — v)j(v) 



due'(u)j(u) + — e(ir) 



Nh 
An 



N h 



It 



2tt 



a 



r 
2 



A/, 

E 



e a 



A,, 



du£(u)j(u) + Ep' 



VP 



(3.8) 
(3.9) 
(3.10) 



Equation ( |3.8| ) can be solved by Fourier transform, paying attention to the fact that j(u) is 



not periodic, but obeys the quasiperiodicity condition j(w + 27r) —j(u) = — Alternatively. 



and with identical results, the symmetric integral operator (1 + 2~K) at the left side of ( |3.8| ) 
can be formally inverted and the solution plugged in ( |3.9| ) and ( |3.10| ) pT| . The result has the 
usual additive form 



A';, 



AE = 2sinh7Ee( 



a 



(3.11) 



where the dressed energy e{u) satisfies 
1 r n 

e(ti) + — — / dvK[u — v)e{v ) = £(w) 
and therefore coincides with R{u) = ^-dn \ k ). As to the momentum, one has 



A,, 



p(af^) + ih 



(3.12) 



where p(a) has been defined in Q2.17 ). The calculation for r odd differs slightly in the 
intermediate steps but also yields ( ^.11 ) and ( ^.12 ), which are then true for Nh arbitrary 
(but obviously even). The fact that Nh is even was missed in [B^] where, following the same 
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assumption made in [0|, one hole was kept fixed at the edge. In other words, only a subset 
of the 2-hole band of states was dealt with 0. 

Eq. (|3.11|) and ( 3.12 ) are simple generalizations of their limit at h = 0, since h appears 



only additively in AP and, implicitly, in the position of the hole which is bound to be on 
the curve C. This simple dependence could not be derived immediately from ( [2.5|) and ( [2.7|) 
because the rapidities {a} depend on h in a nontrivial way through (|2.3|) , and only the explicit 
calculation guarantees that (|3.11|) and Q3.12|) are correct. 



Several comments are in order. Unlike the energy, which being the eigenvalue of a non- 
hermitian operator can, and indeed does have an imaginary part, the momentum must be 
real. This is guaranteed by ( |2.19| ), since E C. It is well known that, from ( |2.7|) , ( |2.8| ) and 



the oddness of 0(a), the momentum can be obtained by summing (|2.8| ) over k 



c\ n 

P = -JfE I » ( 3 - 13 ) 

k=l 

The momentum of the ground state P° is therefore always zero, while the momentum of an 
excited state is 



2tt Nh 



(h) 

N 



k=l 

from which one sees that, as iV — > oo, — | < AP(hole) < |. This is also confirmed by ( [2.171) 



and (|3.12|) . The dispersion relations are obtained by eliminating in (|3.11|) and ( p,12|) 
using 



dn(a; k) = y 1 — k 2 sin (am(a; k) 
which yields 



AE(AP) = m 0] jl-k 2 sin 2 (AP - ih^j m = 2sinh7^^ (3.14) 

An apparent discrepancy with Gaudin's result (for the case h — 0) is clarified in appendix 
B. The dn(a) function (and consequently e(a)) has a non negative real part in the rectangle 
2 One of the authors (GA) is grateful to prof. C. Destri for pointing this out. 
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[—it; — 7r — «7; ir — 27; tt] and this lifts the ambiguity in the sign of ( |3.14| ). It also confirms 
that the choice of the ground state was correct, because the real part of the energy, at least 
under 'small' variations (a countable number of holes), increases. The mininum of Re(AE) 
is reached when = A or B. When this happens AP = ±| and the gap in the spectrum 
is (remember that 2 holes are present in the lowest excited state) 



AP(gap) = 2m \Jl - k 2 cosh 2 (/i) (3.15) 

which guarantees that the 'mass gap' is real (see table 3 for a comparison with numerical 
results). In particular it vanishes at b = —7, that is when A = —71 — ij,B — ir — ij and 

ft = ftc = 2 + E (- r sinh ^> 



2 ^0 ^ cosh (727) 
This is most easily seen from ( |3.11|) 



AE(gap) = 2m dn (^-(±tx - 17); k^j = 2m dn(±I(k) - il'(k); k) = 

Therefore, from (|3.15|) an alternative equation for h c is 

cosh(/i c ) = — 
k 

It is particularly interesting to see how the mass gap vanishes as h — > h c 

AE(gap) ~ 2lm VF(/i c - h)* + 0(h c - h) (3.16) 

The vanishing with an exponent ~ is peculiar of the point under consideration, as it will 
appear clear from the general case, to be discussed later, which includes the vertical field V . 
Finally, we specialize ( |3.15| ) at h — h c . Then the hole excitations are massless and if we set 



AP = — I + e or AP = | — e, < e << 1, we get, respectively 

AP(AP) ~ 2m V=F2ifc' (3.17) 

This dispersion relation is certainly surprising and reflects itself in the peculiar behavior of 
the finite size corrections of the low-lying energy gaps at h — h c . In marked contrast with 
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the 0(jj) scaling typical of spin chains which describe conformally invariant models in the 
continuum limit P3| , we find 



AE ~ + OIN- 1 ) 

where c depends on the state under consideration. The momentum being quantized in units 
of t£ on the finite lattice ( |3.13| ), this behavior is well in agreement with (|3.17j) . 



The sector S z = deserves a special comment. Not knowing what takes the place of 
strings, the analysis of the excitations has been necessarily numerical. Two things have been 
determined. Setting E (S Z = 0,N,h) and Ei(S z = 0,N,h) to be respectively the ground 
state and the lowest-lying of the first band of excited states in the sector S z = 0, on a chain 
of N sites, and at fixed horizontal field h < h c , we found 



AE(S Z = 0,h)= lim 



E X {S Z = 0, N, h) - E (S Z = 0, N, h) 
to be positive, non-zero for h < h c and 



lim AE(S Z = 0,h) = 

h—*h c 

so that, even in this sector, the spectrum becomes massless at h — h c (see table 4). Secondly, 
the 0{N~^) scaling is preserved at h c 

Ex{S z = 0, N, h c ) - E (S Z = 0, N, h c ) ~ iV > 1 

Ei must not be confused with the other (degenerate in the thermodynamic limit) ground 
state that appears in this sector at momentum P = n and is responsible for the spontaneous 
breaking of the arrow-reversal symmetry in the symmetric 6-vertex model |8], [T?J. See table 5 
for a comparison with numerical results. Since we do not study the order parameter (staggered 
polarization) this state will not be discussed here. 

We can now reintroduce V, whose effect is to shift the spectra at S z ^ 0. The mass gap 
for a state with n = n ±r, and consequently 2r holes is easily read from ( |2.5|) and (|3.15|) 



A£(gap; n) = 2rm ^l - k 2 coslr (/i) - 2V(±r) (3.18) 
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An alternative way to reach the boundary with the massless phase is to have a sufficiently 
large |V|. From (|3.18|) the crossing occurs at 



V = ±m \Jl - k 2 coslT 0) (3.19) 

and moves the ground state to sectors of S z > (i.e. n < n ) if V > 0, and to sectors of 
S z < (i.e. n > uq) if V < 0, as it was intuitively predictable from ( |2.5|) . Notice that, unlike 
what happens in (|3.16|) , the mass gap goes to zero linearly in V or linearly in h if V ^ were 
kept fixed and h — > h(V), where h(V) is defined by fl3.19| ). The point V — 0, h — h c (or 
equivalently h = —h c ), where the exponent | of ( |3 . 1 6| ) appears, is clearly special. Even if it 
were approached by changing h and V simultaneously, the term (h c — h)* would dominate 
over the linear term in V. There is no way to erase this effect because it is impossible to 
reach (h c , V = 0) by changing V only: the line h = h c in the (h, V) plane is tangent to the 
phase boundary curve defined by ( |3.19|) . 

This result may look odd, because the energy difference between sectors of different S z 
corresponds to the step free energy for the statistical model, and from (|2.4|) it is hard to see 
how it could vanish other than linearly. Yet it is readily seen that the phenomenon is not an 
artifact of the spin chain. An explicit calculation of the step free energy 



step 



lnA max (5 z = l)-lnA max (5 z = 0) 



can be bypassed observing that the vanishing of / step signals the transition to the incommen- 
surate phase and therefore must be given by ( |2.25| ) 

/step = 2w - 7 + |7 - 2u - b\ 



- (-r sinhK 7 -l T -2^-, m _ 7 < fe < 7 
n coshrry 

The points (h c ,v c ) and (—h c , —v c ), reached on T when b = —7 (or 7) are the equivalent of 
(±/i c , V = 0) in the spin chain phase diagram. They, again, cannot be approached from the 
flat phase by changing v only, since the line h = h c in the (h, v) plane is tangent to T. But, 
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from ( |2.18| ), near b = —7 



h r — h 



2 V 



7T 



k\b + 1 y 



and since / ste p is linear in b near b = —7 (unless u=0) 

/step ~ const(h c — h)i + const(v — v c ) 

at (h c ,v c ). This shows that, like for the spin chain, the exponent | dominates and signals 
that the points (h c ,v c ) and (— h c , —v c ) are essentially different from the other points of T. 

As to the sector S z = 0, we have to extend the numerical analysis carried out for the spin 
chain. If A (S Z = 0,N,h,v c ) and Ai(S z = 0,N,h,v c ) are the largest and next-to-largest 
eigenvalues on the finite lattice in the sector under consideration, we find that 



AA{S Z = 0,h,v c 



lim 



A^S' = 0, N, h, v e ) - A (S Z = 0, N, h, v c ) 



is positive for h < h c and vanishes when h = h c . Furthermore 
- InA^S' = 0,N,h c ,v c ) -\nA (S z = 0,N,h c ,v c 



N2 



iV> 1 



in perfect correspondence with the spin chain scaling of low-lying excitations (see table 6). 



4 The exponent § of the free energy singularity 



As the field crosses the critical value of the T line ( |2.18| ), ( |2.25| ) the system enters a phase where 
horizontal and vertical polarizations change continuously This is an incommensurate phase 



belonging to the universality class of the gaussian model [10|. It is interesting to determine 
the singularity of the free energy as (h, v) approach V from the incommensurate regime. It is 
widely believed [§, |§ that the free energy singularity should be governed by an exponent |, 
but an exact calculation has been done by Lieb and Wu when h = only ||, in which case 



/ ~ c(7,w) 



v - v c (j,u,b = 0) 
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Our calculation is an extension of Lieb's and Wu's method. We will apply it first to the 
ground state energy of the spin chain, and later extend if to the free energy of the statistical 
model. 

Eqs. ( p.9|) and fl2.12|) -( p.l3| ) determine, through the solution of ( [2.11|) , e , y and h as 
functions of A and B. We suppose that such dependence is analytic and e , y and h can be 
expanded in powers of 8 A, 8B as A — > A + 5A, B — > B + 5B. Making explicit the dependence 
of R(a) on A, B by writing R(a; A, B) we have, from ([2. 12 ) 

d A y(A,B) = ~ [ B dad A R(a;A,B) + ^R(A;A,B) 

Z7T J A ZTT 

d B y(A,B) = ~ f B dad B R(a;A,B)-^R(B;A,B) 

Z7T J A Z7T 

5y = d A y{A,B)5A + d B y{A,B)5B + 0{5A 2 ,5B 2 ,5A5B) (4.1) 
Likewise, the energy per site 

1 

e Q (A, B) = cosh 7 + — / dae(a)R(a; A, B) — Vy = cosh 7 — Vy + ejp (A, B) 

Z7l J A 

yields the derivatives 

d A e { o ] (A,B) = — [ B dae(a)d A R(a; A, B) - —e(A)R(A; A, B) 
2n J a 2tt 

1 r B 1 

dae(a)d B R{a; A, B) + —e{B)R(B; A, B) 



d B e£\A,B) 
Se (A,B) 



2-kJa v ' v 7 2tt 



(4.2) 



etc. Similar equations can be obtained for h, specializing ( |2.9| ) to the endpoints of the curve 
and taking the symmetric form 

1 



-4ih(A, B) = p°(A) + p°(B) - — [ d(3R(fc A, B) 

Z7T J A 



Q(A - (3) + Q(B - (3) 



hence 



Aid A h{A,B) = -— ( B d[3d A R([3; A, B) 
2n J a 



®{A-f3) + ®{B-j3) 



+R(A;A,B) 

1 rB 



1 + -e(B-A) 

Zn 



-4id B h(A,B) = -— f d/3d B R((3;A,B) 
Ztt J a 



Q{A- f3) + Q{B - 13) 



+R{B;A,B) 



2tt 



s(b-a) 
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etc.. Equations for the derivatives of R(a; A, B) are readily obtained from 

d A R(a;A,B) + — [ B d(3K(a- (3)d A R(f3;A,B) = —K(a — A)R(A; A, B) 
2tx J a 2tx 

d B R(a; A,B) + ^- [* d$K{a - (3)d B R(P; A, B) = -^-K(a - B)R(B; A, B) 

2,71 J A 2ix 



etc.. 

We have carried out these expansions to the third order in 5 A, 5B. In principle they can be 
used for any A, B with a = ir, \ b\ < 7, and the integrals computed by Fourier transform. How- 
ever, as it is already evident from the first order terms, the expansions simplify considerably 
when carried out around A = —n ± 27, B = n ± 27 which are zeros of the dn function 

R(A ; A , B ) = R(B ; A , B ) = 

The details of the expansion are lengthy but straightforward, so only the final form is of 
interest. Writing 

5A = -5a + i5b 5B = 5a + i5b 

and considering first the expansion around A Q = —n — 27, B = n — 27, we have 

5e = 2c 2 [(5a) 3 -35a(5b) 2 ]-V5y + ... (4.3) 
Sy = -2l5a5b + -5b(Sa) 2 + . . . (4.4) 

7T 7T 

Sh = |[(5a) 2 - {6b)*] + |(5a) 3 + • • • (4-5) 



with 



Cl = A;'(^) 2 >0 

7T 

C2 sinh 7 



1 ^ (— l) n nexp(— 717) 



4 ^ cosh 7 



> 



ci 3tt 

C3 = 1 y^ exp(-|w|7) ^ Q 
ci 7r „ 2 cosh 722 

Notice that 5y = = 5e when 8a = 0, as it should, since by taking 5a = and 5b > we 
reenter the 'flat phase'. It is also important to check that if 5a = there is no way to increase 
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h by changing b, as already discussed in section 2. An increase in h, when keeping y fixed 
at y = 0, can instead be achieved by 5 a < 0, 5b = 0, which confirms the numerical findings 
presented in section 2. A variation 5a > is ruled out a priori, because the periodicity of 
(|2.3| ) in the real direction implies that rapidities are contained in the strip —tt < Re(a) < n 
and a cannot exceed it. 

Another point to discuss is the reliability of (|4.3| )- (|4.5| ) when n > N/2, that is 5y < 0. 
The Bethe-ansatz equations for the symmetric six vertex model are always discussed keeping 
n < N/2, since the Z 2 symmetry of arrow reversal guarantees that the spectrum is the same 
when N/2 < n < N . It is not immediately clear what happens to (|2.3|) when n > N/2. As 
an example, consider the one dimensional sector S z = —N (n = N), whose only eigenstate is 



| Hi ■••!)■ It is not obvious that (|2.3| ) should have only one solution when the number of 
unknowns is N. To be on the safe side we will trust ( |4.3|) -( |4~5l) only for n < y (y > 0). In 
this case 5b > and 5a < 0. To deal with the states at n > y (y < 0), one must resort to 
(|5.4j) which implies 

eo(7> h c + Sh, V, -y) = e (7, -h c - 5h, -V, y) (4.6) 

Hence it is necessary to consider also an expansion around h = —h c ,y = 0. This can be done 
evaluating (fO) and the following equations at the endpoints A' = —n + ij, B' = n + ij. The 



result is that (fO|)-(fO|) still hold, with 5b > (< 0) if 5y > (< 0). A final observation 



about (I~3")-( fO)D is that it is legitimate to neglect higher order terms in ( [4.4| ) and ( ^4.5| ). The 



parameters 5a and 5b are independent and there is no control over their relative magnitude, 
but the second order term in (|4.5|) is dominant unless 5a ~ ±5b, in which case the third order 
term is certainly larger than all possible fourth order terms. Likewise, in ( |4.4| ), no term 5a n 
or 5b n is allowed since we know that 5y = if 5a = or 5b = 0. Consequently, all higher 
order terms can certainly be neglected and one can further limit the expansion to 

5y = — -5a5b 

7T 
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Suppose now that h is kept fixed at h c and V ^ 0. Then, from (p~5*|). 
5b = =F<5a 

where the upper (lower) sign holds for Sy > (< 0). Consider first 5y > 0. Then 
by = — da 

7T 

8e (Sa) = — 4c 2 5a 3 -8a 2 

7T 

which has a minimum, when V > 0, at 

OOo = 

6-7TC2 

that yields 

2 /Vci^ 3 



<yeo(ft = /i c ,y>o) = — 21 ( . 

C2 V 07T 

Notice that no minimum occurs if V < 0. Instead, if we consider < 0, one has a minimum 

at 

8a = 

07TC 2 

that yields 

5e (h = h c ,V<0) = ^(^-) 3 

C 2 V 07T / 

when V < 0. Consequently, 

e (/i = /i C) l/) = e (/i = /i c ,y = 0)-^(|^)V| 3 



is the ground state energy singularity as one approaches the point (h c , V = 0) along the V 
direction. 
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The case V = 0, 8h ^ is more involved. We want 8h > 0, in order to move into the 
incommensurate phase. From 



±Jf(Sa) f(Sa) = Sa 2 + --Sa 3 - —Sh (4.7) 
v 3 c\ c\ 



ci 



7T 



Sy = -^5a(±Jf(5a)) (4.* 



oe (oa) = 2c 2 5a( -2(oa) 2 - 2-(<5a) 3 + ) (4.9) 

Ci Ci 



where the sign in (|4.8| ) depends on whether we want Sy > or 8y < 0. The variation Sa must 
be negative and contained in a range where f(5a) is non negative, so if 

f(5a ) = oa = ~]j^( 5h ) 1/2 + °( Sh ) 
we consider 
5a < Sao 

It is not difficult to see that there is a left neighborhood of Sao (of the order ^/i 1 / 2 ) where 

1. f(5a) is positive 

2. Sy(Sa) is monotonic 

3. o"eo(5a) is decreasing 

Consequently, regardless of the sign in (|4.8|) , Sa = Sa is a local minimum of 5eo(8a), 



which, incidentally, corresponds to Sy = 0. Inserting Sao in ( |4-9|) 



/ 2 \ 3 / 2 

e (/i c + <Wi, V = 0) = e (/i c , V = 0) - 2c 2 ^— <%J (4.10) 

Although it is not obvious from the previous proof, Sy = is actually a stationary point for 
5eo(8y). In fact 



dSy &l 



dSa 'Sao 
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and 

dSy ( ci rjTz— c x f'(Sa) 



dSa 



becomes infinite at Sao. 

The calculation of the free energy singularity is a simple extension of this method. The 
variation of the ground state energy is now replaced by (see (|2.23|) ) 

-S(F(u, 7, h, y) + vy) = -SF - Sy(v c + Sv) 

The variation of F is computed by means of an expansion analogous to (|4.1|) - (|4.2|) . To keep 
things simple we consider d < 7 — 2u, which is certainly true for u sufficiently small, so that 
Aft dominates over A^. The surprisingly simple result is that, like for the spin chain, the first 
nonzero contribution comes at the third order. The quantity to minimize is 

2c' 2 (Sa 3 - 3SaSb 2 ) - SySv 

where 

, _ ^ 1 (-l)"cosh( 7 (n-2n)) 1 
C2 -6n [ 2 + ^ cosln» J>U 



which looks exactly like (|4.3|) - fl4~5|) provided c 2 — > c' 2 and V —> Sv. The conclusions are 
therefore the same and the free energy leading singularities approaching (h c , v c ) from the 
incommensurate phase are 

/ 2 \ 3 /2 

f(u,<y,h c + 5h,v c ) = f(u,<y,h Cl v c )-2c 2 \^—Shj 

2 / Cl ^ 3 



2 ( c\ \ A 

f(u,j,h c ,v c + Sv) = f(u,j,h c ,v c ) - -jzl^—j \Sv 



5 Discussion 

It is interesting to speculate about the nature of the phase transition at (h c , v c ) and compare 
it with what happens at the other points of T. 
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We have always worked in the assumption that a complete set of eigenstates exists for the 
transfer matrix. It is well known then that the correlation function of 2 vertical arrows can 



be analysed through a spectral decomposition [24|. For the correlator of two vertical arrows 



along the same column one has, on a iV x M lattice 



n ^.zrpM — n\ 



(ok|*> 



2 / A , \ n 



(5.1) 



and for the correlation function along the horizontal direction 

Ko «„, ) = Tr( ^g M) ^ (0K<|0) (5.2) 

where we have denoted with |0) the eigenstate of the largest eigenvalue of the transfer matrix 
on a finite lattice of width N. Here Y^,k denotes the sum over a complete set of eigenvectors 
of the transfer matrix. It is useful to consider first what happens at h = 0. Since, obviously 



N 



3=1 







and the ground state lies in the sector S z = 0, only this sector contributes to ( |5.1|) and (|5.2j). 
Furthermore, v does not change the spectrum within a sector of fixed S z , and it does not 
modify the eigenvectors which depend (as explicitly seen from the Bethe Ansatz || [7J]) on 7 
and h only. We conclude that v has no effect whatsoever on the correlators, which remain 
those of the symmetric six-vertex model, until the level-crossing transition takes place at 
v = v(j,u,b = 0) (as given by (|2.25f) ). Here the system moves into the gaussian phase, and 
the ground state even for small 5v = v — v(j, u, b = 0) falls into a sector at y ^ 0, hence 
with S z of order N. The transition occurs without divergence of the correlation length which 
jumps from the (finite) value of the symmetric six-vertex to infinity. 

It is our conjecture that this picture does not change when — h c < h < h c , and the level 
crossing transition persists along all T, up to the points (h c ,v c ) or (— h c , — v c ). Here the 
correlation length should diverge according to the following argument. The fact that in the 
sector S z = 0, as found numerically in section 3, 

lim_ AlnA(S 2 = 0,h,v c ) = 

h — >h c 
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is, by itself, not sufficient to prove the divergence of the correlation length £1 in the vertical 
direction. The fact that in ( |5.1|) are generally complex forces one to sum over a whole band 
of them because oscillations can affect the behavior of £1 f24 |. This in not doable until it is 
clarified what takes the place of strings in the sector under consideration. Nevertheless, one 
can look at the horizontal correlation function in fl5.2|) . As shown in section 4, it is possible 
to enter the gaussian phase at (h c , v c ) keeping the ground state at y = 0. From fl2.10|) y = 



does not necessarily imply that S z = (n = y), since S z can be nonzero but remain finite in 
the limit N — > oo and y would still be zero. Still it is tempting to conjecture that the ground 
state remains at S z = 0, and this is confirmed by preliminary numerical results on the spin 
chain. Hence no level crossing occurs here, because the Perron-Frobenius theorem prevents 
if from happening in a sector of fixed S z , and the expectation value in ( |5.2j ) is taken with the 
same Bethe-Ansatz eigenvector in the 2 phases. Since in the gaussian phase the horizontal 
correlation lenght £2 is infinite, ^(h) must diverge as h — > ±/i c . 

Although these conclusions are rather speculative, they seem to warrant a further investi- 
gation of the phase transition at issue. It should be recalled that, if indeed (h c ,v c ) is critical 
in the sense of diverging correlation length, it cannot be classified in terms of conformal field 
theory as shown by the anomalous scaling discussed in section 3. 
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Appendix A 

We discuss here the symmetries of TC and of the transfer matrix. Under the action of the 
(unitary) charge conjugation operator 

jv 

C = J] a x k C = C f = C- 1 

k=l 

H. transforms as 



CH{j, h, V)C = ft(7, -h, -V) 



(5.3) 



Since 



A? 



N 



3=1 3=1 

the spectra of the sectors at fixed S z are related by 
S{ 1 ,h,V,S z )=S{ 1 ,-h 1 -V,-S z ) 



(5.4) 



The same symmetry operation can be applied to the transfer matrix, using the matrix form 
of C 



N 



{a, Co!) = Y[5, 



k=l 



From ( p.l| ) and the definition of the Botzmann weights it is elementary to see that 



C T(-f,u,h,v) C 



T(7,m, -h, -v) 



a.a_ 



(5.5) 



which implies the relation between spectra 



<S T M{h, v, S z } = SrM{-h, -v, -S z } 



(5.6) 



This symmetry manifests itself in the fact that the partition function 



Z PF (h, v) = Z PF (-h, -v) 
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As far as 7i is concerned, another symmetry is in effect, implemented by the space inversion 
operator (not to be confused with the momentum) 

P\cni, a 2 , ■ ■ ■ , &N-1, ®n) = \®n, ajv-ij • • • j a 2, «i) P 2 = 1 

PH(y,h,V)P = H(7,-h,V) 

Hence 

C P H(j, h, V = 0) P C = H(y, h,V = Q) 

So at V = the spin chain recovers the Z 2 symmetry under spin reversal and the spectra at 
S z and — S z are identical. It is noteworthy that the same is not true for the transfer matrix. 
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Our Hamiltonian, at h = 0, V = 0, can be written, neglecting a constant additive term 

1 N 



i=i 



A ^Z—Z i i _J/ _J/ 



with A < 1 



(5.7) 



and it is mapped onto the Hamiltonian (see, e.g. pT| ' 



N 



with A > 1 



(5.f 



through a unitary transformation 

AT/2 

H G = UHU- 1 U=H &2j—i 

i=i 

Consider the shift operator S of ( |2.6| ). Let |\l/ n , P) be an eigenstate of H with momentum P 
and S 12 = y — n 



N 



N 



-E a] |* n ,P) = (--n) |* BJ P> 



5 |* n ,P) = exp(-zP) |* n ,P) 
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Define Uq as 

N 

U |* n ,P):=n°i l^n,P) = (-l) n l^n,P) 

i=i 

Then, an eigenstate of Tia with the same energy is U |\I/„, P). Since 
S U |tf„,P) = exp(-tP) S U S- 1 \V n ,P) 

N/2 

= exp(-iP) n °2i \*n,P) 

exp(-iP) U U \V n ,P) = exp(-i(P + n7r)) U |* n ,P) 

= exp(-zP) [/ |* ft ) 

£/ |\I/ n ,P) has momentum P = P + ror. Hence the ground state of (|5.8| ) has momentum 
P = n 7T = y7r, while P = 0. As to the excitations 

AP = AP + (n - n )n = AP - 

Each hole carries an additional momentum — 71/ 2 which implies 

AP = m yjl- Psin 2 (AP) -> AP = m ^l - fc 2 cos 2 (AP) 
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Lattice size 


a 





8 


2.226901585 


-1.186380535 


12 


2.522418777 


-1.215745628 


16 


2.674453171 


-1.227403536 


20 


2.766810718 


-1.233071106 


24 


2.828778340 


-1.236225697 


28 


2.873206253 


-1.238154458 


32 


2.906605731 


-1.239417236 


36 


2.932624180 


-1.240288031 


40 


2.953462078 


-1.240913442 


Extrap. oo 


3.141592659(9) 


-1.243603723(1) 




7T = 3.141592653 


- 7 = -3.737102242 



Table 1: Position of the endpoints ±a + ib of the curve C for h < h c ; 
exp(2h) = 3 
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Lattice size 


a 


I. 



8 


1.345603261 


-3.0497549071 
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-3.2366686985 


16 


1.535932819 


-3.3472122872 


20 


1.564274277 


-3.4190538006 


24 


1.580405494 


-3.4690826450 


28 


1.590409443 


-3.5057653322 


32 


1.597021761 


-3.5337453681 


36 


1.601612009 


-3.5557588049 


40 


1.604924714 


-3.5735136386 


Extra p. oo 


1.6193362(3) 


-3.737102232(7) 




n = 3.141592653 


- 7 = -3.737102242 



Table 2: Position of the endpoints ±a + ib of the curve C for h > h c ; 
exp(2/i) = 18 
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Lattice size 


e = l.U 


e = 9. (J 


8 


38.8549946261240 


27.3592680025296 


12 


38.4729554492734 


23.5353424072427 


16 


38.3095545699746 


21.4216001709441 


20 


38.2251709688714 


20.0937869680560 


24 


38.1760618666305 


19.1948069770300 


28 


38.1450125441243 


18.5549980007490 


32 


38.1241500821837 


18.0828136699350 


36 


38.1094640015976 


17.7244227601235 


40 


38.0987382252780 


17.4461598430727 


Extrap. oo 


38.04991359(1) 


15.8884(8) 


Exact eq. 


38.04991361 


15.8887254 



Table 3: Energy gap of the spin chain for the first excited state for different values of h 
compared to the analytical result (3.15) 



36 



Lattice size 


J1K n n 

e = 9.U 


e = 9.5 


9/i inn 

e = 1U.U 


~2hr in r 1 TOT 

e^" c = 1(J. 51787 


8 


40.28463586426577 


40.60665243638542 


41.00560893816288 


41.50186071099847 


12 


32.50806556257382 


32.33242547913501 


32.27369006467646 


32.34509918439876 


16 


28.03119022066448 


27.47527007849106 


27.06291429820305 


26.80804074773090 


20 


25.23278497081201 


24.37422080701744 


23.67516261790179 


23.15771126310622 


24 


23.35422673956813 
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20.56522041772249 


28 


22.02474072486461 


20.70732422258583 


19.55338243467817 


18.62033939778466 


32 


21.04552971555199 


19.54981005828927 


18.20747202515815 


17.10061549106331 


36 


20.30141022559842 


18.65287555243093 


17.14104807658379 


15.87553491911761 


40 


19.72145517781357 


17.94130164070793 


16.27579201027032 


14.86352304583235 


Extra p. oo 


15.946(3) 


13.054(9) 


9.37499(4) 


-1.3(3) • 10~ 5 



Table 4: Mass gap of the spin chain for different values of h < h c 
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Lattice size 


Exponent 


8 


-0.6147896371667922 


12 


-0.6526671750014708 


16 


-0.6559635919155932 


20 


-0.6511929107450443 


24 


-0.6444786397317260 


28 


-0.6376031966136601 


32 


-0.6311198906745727 


36 


-0.6251786621408798 


40 


-0.6197881465618827 


Extrap. oo 


-0.5031(1) 


expected value 


-0.5 



Table 5: Exponent of N for the finite-size corrections of the energy gap in the sector S z 
on the critical line h = h c 
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Lattice size 


Exponent 


8 


-0.6696077334198590 


12 


-0.6362863326795953 


16 


-0.6165741467010993 


20 


-0.603244670912252 


24 


-0.5934852417528157 


28 


-0.5859536198128706 


32 


-0.5799199948938484 


36 


-0.5749497141578801 


40 


-0.5707659045338731 


Extra p. oo 


-0.499998(4) 



Table 6: Exponent for the scaling of the free energy gap in the sector S z 
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Figure 1: Boltzmann weights in the notation with spectral parameter u compared to that of 
ref. [9]. The physical region is < u < 7 
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Figure 2: Distribution of the rapidities a k for iV = 44. The lower curve corresponds to the 



ground state in the sector with n = y and exp(2/i) = 9, the upper curve to the ground state 
of the sector with n = 4 and exp(2/i) = | 
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Figure 4: Distribution of I^s for the ground state and an excited state with 
rapidities 
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